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ABSTRACT 

We are using the Very Long Baseline Array and the Japanese VLBI Ex- 
ploration of Radio Astronomy project to measure trigonometric parallaxes and 
proper motions of masers found in high-mass star-forming regions across the 
Milky Way. Early results from 18 sources locate several spiral arms. The Perseus 
spiral arm has a pitch angle of 16° ± 3°, which favors four rather than two spi- 
ral arms for the Galaxy. Combining positions, distances, proper motions, and 
radial velocities yields complete 3-dimensional kinematic information. We find 
that star forming regions on average are orbiting the Galaxy ^ 15 km s^^ slower 
than expected for circular orbits. By fitting the measurements to a model of 
the Galaxy, we estimate the distance to the Galactic center Rq = 8.4 ± 0.6 kpc 
and a circular rotation speed Gq = 254 ± 16 km s~^. The ratio O0/-R0 can be 
determined to higher accuracy than either parameter individually, and we find it 
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to be 30.3 ± 0.9 km kpc~^, in good agreement with the angular rotation rate 
determined from the proper motion of Sgr A*. The data favor a rotation curve 
for the Galaxy that is nearly fiat or slightly rising with Galactocentric distance. 
Kinematic distances are generally too large, sometimes by factors greater than 
two; they can be brought into better agreement with the trigonometric parallaxes 
by increasing Qq/Rq from the lAU recommended value of 25.9 km s~^ kpc~^ to 
a value near 30 km s~^ kpc^^. We offer a "revised" prescription for calculating 
kinematic distances and their uncertainties, as well as a new approach for defin- 
ing Galactic coordinates. Finally, our estimates of 6o and GoZ-Ro, when coupled 
with direct estimates of Rq, provide evidence that the rotation curve of the Milky 
Way is similar to that of the Andromeda galaxy, suggesting that the dark matter 
halos of these two dominant Local Group galaxy are comparably massive. 

Subject headings: Galaxy: fundamental parameters, structure, kinematics and 
dynamics, halo — stars: formation — astrometry 



1. Introduction 

The Milky Way is known to possess spiral structure. However, revealing the nature 



of this structure has proven to be elusive for decades. The iGeorgelin fc Georgelin I (1l976l ) 



study of HII regions produced what has been generally considered the "standard model" for 



on this structure. Indeed, we c 


io not really know the number of s 
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) or how 



tightly wound is their pattern. The primary reason for the difficulty is the lack of accurate 
distance measurements throughout the Galaxy. Photometric distances are prone to calibra- 
tion problems, which become especially severe when looking through copious dust to distant 
objects in the plane of the Galaxy. Thus, most attempts to map the Galaxy rely on ra- 
dio frequency observations and kinematic distances, which involve matching source Doppler 
shifts with those expected from a model of Galactic rotation. However, because of distance 
ambiguities in the first and fourth quadrants (where most of the spiral arms are found) and 



non-circular motions, kinematic dist 
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We are measuring trigonometric parallaxes and proper motions of sources of maser 
emission associated with high-mass star forming regions (HMSFRs), using the National Ra- 
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dio Astronomy Observatory's Very Long Baseline Array (VLBA) and the Japanese VLBI 
Exploration of Radio Astronomy (VERA) project. The great advantage of trigonometric 
parallaxes and proper motions is that one determines source distances directly and geomet- 
rically, with no assumptions about luminosity, extinction, metallicity, crowding, etc. Also, 
from the same measurements, one determines proper motions, and if the time sampling is 
optimal there is little if any correlation between the parallax and proper motion estimates. 
Thus, the magnitude of the proper motion does not affect the parallax accuracy. Combining 
all of the observational data yields the full 3-dimensional locations and velocity vectors of 
the sources. 

Results for 12 GHz methanol (CH3OH) masers toward 10 HMSFRs, ca rried out with the 



VLBA (program BRIOO), are reported in the first five papers in this series (IReid et al. 112009 



MoscadeUi et al. I l2009l : IXu et al. I l2009l : IZhang et al. 



20091 : iBrunthaler et al. 1 12009^ . here 



after Papers I through V, respectively. Eight other sources with H2O or SiO masers have been 



measured with VERA (IHonma et al. 1120071 : iHirota et al. 1120071 : IChoi et al. 



2006 



2008 ) and with methanol, H2O or continuum emission with the VLBA (jHachisuka et al. 



Hachisuka et al. 



Menten et al. 2007; Moellenbrock. Claussen fc Goss 2007; Bartkiewicz et al. 



2008 



Sato et al. 



2008 



20091 ). In this paper, we collect these parallaxes and proper motions in 



order to study the spiral structure of the Galaxy. Combining positions, distances, Doppler 
shifts, and proper motions, allows us not only to locate the HMSFRs that harbor the target 
maser sources in 3-dimensions, but also to determine their 3-dimensional space motions. In 
§5] we map the locations of the HMSFRs and measure the pitch angles of some spiral arms. 
In 113] we use the full 3-dimensional spatial and kinematic information to examine the non- 
circular (peculiar) motions of these star forming regions. We also fit the data with a model 
of the Galaxy and estimate the distance from the Sun to the Galactic center (i?o) and the 
circular orbital speed at the Sun (60) • The nature of the rotation curve and its effect on 
estimates of Rq and 60 is also discussed. In ^we compare kinematic distances with those 
determined by trigonometric parallax and offer a new prescription to improve such distance 
estimates. In ^ we discuss limitations of the current definition of Galactic coordinates and 
suggest a new system based partly on dynamical information. Finally, we discuss the broader 
implications of our results in §3 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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2. Galactic Spiral Structure 

Table [1] summarizes the parallax and proper motions of 18 regions of high-mass star 
formation measured with VLBl techniques. The locations of these star forming regions in 
the Galaxy are shown in Figure [H superposed on an artist's conception of the Milky Way. 
Distance errors are indicated with error bars (Icr), but for most sources the error bars are 
smaller than the dots. 
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Table 1. Parallaxes & Proper Motions of High- mass Star Forming Regions 



Source 






b 






Parallax 










fly 




^'LSR 


Ref. 




(de 


g) 


(de 


?) 




(mas) 


(mas y ^) 


(mas y ^) 


(km s~^) 




G 23.0-0.4... 


23 


01 


-0 


41 





218 + 0.017 


— 1 


72 + 


04 


—4 


12 + 


30 


+81 + 3 


V 


G 23.4-0.2... 


23 


44 


-0 


18 





170 + 0.032 


— 1 


93 + 


10 


—4 


11 + 


07 


+97 + 3 


V 


G 23.6-0.1... 


23 


66 


-0 


13 





313 + 0.039 


— 1 


32 + 


02 


—2 


96 + 


03 


+83 + 3 


1 


G 35.2-0.7... 


35 


20 


-0 


74 





456 + 0.045 


— 


18 + 


06 


—3 


63 + 


11 


+28 + 3 


IV 


G 35.2-1.7... 


35 


20 


-1 


74 





306 + 0.045 


— 


71+0 


05 


—3 


61 + 


17 


+42 + 3 


IV 


W 51 IRS 2.. 


49 


49 


-0 


37 





195 + 0.071 


-2 


49 + 


08 


-5 


51+0 


11 


+56 + 3 


III 


G 59.7+0.1... 


59 


78 


+0 


06 





463 + 0.020 


-1 


65 + 


03 


-5 


12 + 


08 


+27 + 3 


III 


Cep A 


109 


87 


+2 


11 


1 


430 + 0.080 


+0 


50 + 1 


10 


-3 


70 + 


20 


-10 + 5 


II 


NGC 7538.... 


111 


54 


+0 


78 





378 + 0.017 


-2 


45 + 


03 


-2 


44 + 


06 


-57 + 3 


II 


IRAS 00420.. 


122 


02 


-7 


07 





470 + 0.020 


-1 


99 + 


07 


-1 


62 + 


05 


-44 + 5 


2 


NGC 281 


123 


07 


-6 


31 





355 + 0.030 


-2 


63 + 


05 


-1 


86 + 


08 


-31 + 5 


3 


W3(0H) 


133 


95 


+ 1 


06 





512 + 0.010 


-1 


20 + 


20 


-0 


15 + 


20 


-45 + 3 


4 


WB 89-437... 


135 


28 


+2 


80 





167 + 0.006 


-1 


27 + 


50 


+0 


82 + 


05 


-72 + 3 


5 


S 252 


188 


95 


+0 


89 





476 + 0.006 


+0 


02 + 


01 


-2 


02 + 


04 


+11 + 3 


I 


S 269 


196 


45 


-1 


68 





189 + 0.016 


-0 


42 + 


02 


-0 


12 + 


08 


+20 + 3 


6 


Orion 


209 


01 


-19 


38 


2 


425 + 0.035 


+3 


30 + 1 


00 


+0 


10 + 1 


00 


+10 + 5 


7 


G 232.6+1.0.. 


232 


62 


+1 


00 





596 + 0.035 


-2 


17 + 


06 


+2 


09 + 


46 


+23 + 3 


I 


VY CMa 


239 


35 


-5 


06 





876 + 0.076 


-3 


24 + 


16 


+2 


06 + 


60 


+ 18 + 3 


8 



Note. — Columns 2 and 3 give Galactic longitude and latitude, respectively. Columns 5 and 6 are proper 
motions in the eastward {^x — HaCOsS) and northward directions (fiy — fis), respectively. Column 7 lists Local 
Standard of Rest velocity components; these can be converted to a heliocentric frame as described in the Ap- 
pendix. References are I: iReid et al. I (|2009l ): II: iMoscadelli et al. I |2009[ ): III: IXu et al. I |2009l ): IV: Izhang et al" 



I2OO9I) : V: iBrunthaler et al. I |2009t) : 1: iBartkiewicz et al. I |2008[ ): 2: iMoellenbrock. Claussen k Goss I |2007h : 3: 
Sato et al. I |2008h : 4: IXu et al. I |2006h : iHachisuka et al. I |2006h : 5: iHachisuka et al. I |2009l) : 6: iHonma et al 



|2007h : 7: iHirota et ahl |2007[ l: iMenten et al~l |2007[ l: 8: Ichoi et al. I |2008h . The calculations in this paper use an 
early parallax and proper motion estimate for IRAS 00420+5530 cited above; the values reported more recently 
by IMoellenbrock. Claussen fc Goss I (|2009( ) are slightly different but would not substantively change the results 
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presented here. 
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2.1. Spiral Arms 



The HMSFRs with parallaxes locate several spiral arms. The three sources closest 
to the Galactic center (G 23.0—0.4, G 23.4—0.2, and G 23.6—0.1) appear to be members 
of the Crux-Scutum or possibly the Norma or the 3-kpc arm. However, the parallax un- 
certainties for these distant, low-declination sources are currently not adequate to clearly 
distinguish among these arms, especially in the crowded region where the Galactic bar (see 
Blitz fc Spergdl Jl991bh and re ferences therein) ends and the arms begin (IBenjamin et al. 



2005 



Dame k Thaddeus l[2008l l 



Three sources (G 35.2—0.7, G 35.2—1.7 and W 51 IRS 2) are probably in the Carina- 
Sagittarius arm, whose distance from the Sun is 2.5 kpc at Galactic longitude, £, of ~ 35°. 

Five sources (S 252, W3(0H), IRAS 00420+5530, NGC 281, and NGC 7538) clearly 
trace a portion of the Perseus arm, which is located between distances of 2.10 kpc ai i = 189° 
(S 252) and 2.64 kpc at £ = 112° (NGC 7538). NGC 281 is slightly offset from the other 
sources in the Per seus arm and is believed to be associated with an expanding super-bubble 
( ISato et al. II2008I ). As such, it may not accurately trace spiral structure. 



Two sources, (S 269 and WB 89-437), measured bv lHonma et al. I (120071 ) with the VERA 
array and iHachisuka et al. I (120091 ) with the VLBA, lie beyond the Perseus arm and begin to 
trace an Outer (Cygnus) arm at a distance from the Sun of 5.3 kpc at £ = 196° for S 269 to 
5.9 kpc at £ = 135° for WB 89-437. 

The remaining five sources (G 232.6+1.0, VY CMa, the Orion Nebula, Cep A and 
G59. 7+0.1) trace the Local (Orion) "arm," which appears to be a spur between the Carina- 
Sagittarius and Perseus arms. The Sun is in or near this spur, and we can trace it between 
G 59.7+0.1 near the Carina-Sagittarius arm at £ = 60° and G 232.6+1.0 near the Perseus 
arm at £ = 233°. 



2.2. Pitch Angles 

Spiral arm pitch angles can be estimated when two or more sources can be confidently 
identified as members of a single arm. The pitch angle, ipy of an arm segment can be defined 
by constructing a line segment between sources in the same section of the arm. Next, 
construct a line tangential to a Galactocentric circle that passes through the midpoint of 
this segment and determine the angle between these two lines. An ideal log-periodic spiral 
arm can be defined by the equation 

In {R/Rref) = -(/9 - Pref) t^^lp , 
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Fig. 1. — Locations of high- mass star forming regions for which trigonometric parallaxes have 
been measured. Parallaxes from 12 GHz methanol masers are indicated with dark blue dots and 
those from H2O and SiO masers or continuum emission (Orion) are indicated with light green 
dots. Distance error bars are indicated, but most are smaller than the dots. The Galactic center 
{red asterisk) is at (0,0) and the Sun {red Sun symbol) at (0,8.5). The background is an artist's 
conception of Milky Way (R. Hurt: NASA/JPL-Caltech/SSC) viewed from the north Galactic pole 
from which the Galaxy rotates clockwise. The artist's image has been scaled to place the HMSFRs 
in the spiral arms, some of which are labeled. 
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where R is the Galactocentric radius at a Galactocentric longitude (3 (defined as toward 
the Sun and increasing with Galactic longitude) for an arm with a reference radius Rref at 
Pref- In Figure m we plot log(-R/l kpc) versus (3 for the three arms where we can clearly 
identify two or more HMSFRs. In such a plot, spiral arm sections appear as straight lines. 
Some of the data deviate from fitted lines by considerably more than the parallax errors, as 
expected for variations of the locations of star forming regions within an arm whose width 
is ~ 100 pc. Thus we used unweighted straight line fits to estimate spiral arm pitch angles. 

The Perseus arm sources (excluding NGC 281) indicate a pitch angle of 16?5±3?1 be- 
tween Galactic longitude 112° and 189°. This is in t he uppe r half o f pitch angle estimates of 



5° to 21° for spiral arms in the Galaxy collected by I Vallee I (119951 ). The five HMSFRs with 
parallaxes that trace the Local "arm" indicate a mean pitch angle of 27?8±4?7. However, 
the Local arm is probably not a global spiral arm; instead it appears to be a short segment 
or spur between the Carina-Sagittarius and Perseus arms. 

Two sources (S 269 & WB89-437) appear to be part of the Outer arm and formally 
yield a pitch angle of 2?3. This suggests that the Outer arm might have a smaller pitch 
angle than the Perseus arm. This may be of significance, but with only two sources, more 
parallaxes are needed before reaching any firm conclusions. 

For other spiral arms, we have too few parallaxes to reliably determine pitch angles. The 
sources that are possible members of the Carina-Sagittarius arm (G 34.2—0.7, G 34.2—1.7 
& W 51 IRS 2) would formally give a wide range of pitch angles. However, because one or 
more sources might be associated with the Local arm, we cannot reliably estimate the pitch 
angle of the Carina-Sagittarius arm at this time. 



3. Galactic Dynamics 

Given measurements of position, parallax, proper motion and Doppler shift, one has 
complete three-dimensional location and velocity vectors relative to the Sun. One can then 
construct a model of the Milky Way and adjust the model parameters to best match the 
data. We model the Milky Way as a disk rotating with speed Q{R) = ©o + ^ (-R ~ -^o); 
where Rq is the distance from the Sun to the Galactic center. We started the fitting process 
by assuming a fiat rotation curve (i.e. ^ = 0). Later, we relaxed this assumption and solved 
for followed by an investigation of other forms of the Galactic rotation curve. Since all 
measured motions are relative to the Sun, we need to remove the peculiar (non-circular) 
motion of the Sun, which is parameterized by Uq toward the Galactic center, Vq in the 
direction of Galactic rotation, and Wq towards the north Galactic pole (NGP). Table [2] 
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1.2 




Fig. 2. — Spiral arm pitch angles. The logarithm of Galactocentric radius, R, (in kpc units) is 
plotted against Galactocentric longitude Data based on trigonometric parallaxes for sources 
that can be confidently assigned to a spiral arm are shown along with la uncertainties. Positional 
variations of star forming regions within an arm are clearly greater than the parallax uncertainties. 
Unweighted fits of straight lines to the data are shown with solid lines. Pitch angles are proportional 
to the negative of the arctangent of the line slopes. 
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summarizes these and other parameters. 

We adjusted the Galactic parameters so as to best match the data to the spatial- 
kinematic model in a least-squares sense. For each source, we treated the measured parallax 
(vTs), two components of proper motion {fi^, /iy), and the hehocentric velocity ivHeiio) as data. 
The observed source coordinates are known to extremely high accuracy and were treated as 
independent variables. The model is a smoothly rotating galaxy given by the parameters 
listed in Table [21 Specifically, the model parallax is calculated from a kinematic distance, 
based on the observed Doppler shift. The three-dimensions of motion relative to the Sun 
(proper motion and heliocentric Doppler shift) are calculated from the source location, taking 
into account the size and rotation curve of the galaxy model, and th e solar and source peculiar 



motions. We adopt the Hipparcos determination of solar motion (iDehnen &: Binney Ill998l ) 
as definitive and generally did not vary these parameters. However, in one least-squares fit, 
we solved for these parameters for illustrative purposes in order to compare Solar Motion 
results from Hipparcos stars and our HMSFRs. 

Our choice of weights for the data in the least-squares fitting process requires some 
comment. While the heliocentric velocity of any maser spot can be measured with very 
high accuracy, it may not exactly refiect the motion of the HMSFR. The i nternal motions of 



meth anol masers are generally small and cause uncertainty of ~ 3 km s ^ (IMoscadelli et al. 



20021 ). whereas H2O masers can be associated with fast outfiow and, if not accurately mod- 
eled, can lead to larger uncertainty in the motion of the exciting star. In addition, the Virial 
motion of an individual massive star (associated with the masers) with respect to the entire 
HMSFR is likely to be 7 km s~^ per coordinate (e.g. for a region of mass of ~ 3 x 10^ M© 
and radius of ~ 1 pc). Therefore, we allow for a deviation of the measured motion from the 
center of mass of its associated HMSFR by adding an uncertainty of avir = 7 km s~^ in 
quadrature with the internal motion estimates (between 3 and 5 km s~^). Specifically, the 
variance weights for the wlsr data, w(fLSR), are calculated from w(fLSR) = V(^^lsr ~^'^vir)- 

Since the parallax data is compared to a kinematic model, we considered both the paral- 
lax measurement uncertainty and a modeling uncertainty for the kinematic distance, akd, ow- 
ing to the total uncertainty in the heliocentric velocity of the associated HMSFR. These two 
components were added in quadrature when calculating the weights: w^Hs) = l/(cr^+c"L/^s)' 
where dg = l/vr^. Similarly, the proper motion weights allowed for both measurement uncer- 
tainties and the possible deviation of the measured maser motions from the center of mass of 
the HMSFR. The latter term was set by the uncertainty in the heliocentric velocity divided 
by the distance. Thus, for either proper motion component w(/i) = l/(o"^ + CTy^^/dl). 



Table 2. Galaxy Model Parameter Definitions 



Parameter Definition 

i?o Distance of Sun from GC 

9o Rotation Speed of Galaxy at Rq 

g Derivative of 6 with R: e{R) = Qq + ^ {R - Ro) 

Uq Solar motion toward GC 

Vq Solar motion in direction of Galactic rotation 

Wq Solar motion toward NGP 

Us Average source peculiar motion toward GC 

Vs Average source peculiar motion in direction of Galactic rotation 

Ws Average source peculiar motion toward NGP 



Note. — GC is the Galactic Center and NGP is the North Galactic Pole. The average 
source peculiar motions {Us,Vs,Ws) are defined at the location of the source and are rotated 
with respect to the solar motion {Uq,Vq,Wq) by the Galactocentric longitude, /?, of the source 
(see Appendix Figured]). We solve for the magnitude of each component, but the orientation 
of the vector for each source depends on location in the Galaxy. 
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3.1. Galactic 3-D Motions 



We first used all 18 sources listed in Tabled] and solved only for the fundamental Galactic 
parameters, yielding Rq = 8.2 kpc and 0o = 265 km s~^ for a fiat rotation curve = 0) 
(see Fit 1 in Table [3]). The value of 263 for 70 degrees of freedom was quite large, 
and the post-fit residuals showed clear systematic deviations, indicating a deficiency in this 
2-parameter model for Galactic dynamics. 

Figure [3] shows the peculiar motions of the HMSFRs in the Galactic plane by trans- 
forming to a reference frame that rotates with the Galaxy. Peculiar motions relative to 
two Galactic rotation models are shown, one for Oq = 220 km s~^ (the lAU recommended 
value) and the other for Oq = 254 km s^^ (our best fit value from Section [3^ below) . Both 
transformations assume Rqj= 8.5 kpc, a fiat rotation curve, and the Hipparcos solar motion 



of iDehnen &: Binney I (119981 ). (The equations used for the transformation are documented 
in the Appendix.) Sizeable systematic motions are clearly evident — almost all sources have 
a significant component of peculiar motion counter to Galactic rotation. On average these 
star forming regions orbit the Galaxy 15 km s~^ slower than the Galaxy spins. As is 
evident from the two sets of peculiar motions, this conclusion is insensiti ve to the values 
adopted for Gq. Similarly, adopting a more complex rotation curve, e.g. the IClemens I ( 119851 ) 
curve, would not change the qualitative nature of the residuals. HMSFRs appear to orbit the 
Galaxy slower than for circular orbits. This might be explained by star formation triggered 
by the encounter of molecular gas with a shock front associated with a trailing spiral a rm 
and may help explain the 17 km s~^ dispersion seen in HI data by iBrand &: Blitz I (119931 ). 



For the distribution of sources in our sample, the solar motion parameters Uq and Vq can 
partially mimic the average source peculiar motions. W e belie ve the solar motion parameters 
determined from Hipparcos data bv lDehnen fc Binney I (119981) are well determined, and they 
have been independently confirmed by lMendez et al. I (119991 ). based on the Southern Proper- 
Motion program data. However, it is instructive to solve for the solar motion parameters with 
the parallax and proper motion data. Doing so we find an acceptable fit with Rq = 8.4 kpc, 
Qo = 242 km s-\ Uq = 9 km s-\ Vq = 20 km s"^ and Wq = 10 km s-^ (The value for 
this fit was 67.2 for 59 degrees of freedom, which is somewhat worse than the value of 65.7 
found in §3.21 where we adopt the Hipparcos solar motion parameters and solved instead for 
average source peculiar motion components.) 



In FigureHlwe reproduce the Hipparcos solar motion data from Figure 4 of lDehnen fc Binney 



( 119981 ). Their data were binned by stellar colors, plotted against stellar dispersion and the 
solar motion components estimated as minus the average velocity of all stars in each bin. We 
have also plotted our estimates of the solar motion, plotted at near-zero "stellar dispersion" 
appropriate for newly formed stars. Also included in the bottom panel is the value of Wq 
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Fig. 3. — Peculiar motion vectors of high mass star forming regions (superposed on an artist 
conception) projected on the Galactic plane after transforming to a reference frame rotating with 
the Galaxy. A 10 km s^^ motion scale is in the lower left. The Galaxy is viewed from the north 
Galactic pole and rotates clockwise. The light (yellow) arrows are for lAU standard values of 
i?o = 8.5 kpc and ©o = 220 km s~^ and a flat rotation curve, whereas black arrows are for 
00 = 254 km s~^. This demonstrates that the qualitative result that high mass star forming 
regions orbit the Galaxy slower than the Galaxy rotates is not sensitive to the value of Qq. 
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determined from the proper motion of Sgr A* by iReid fc Brunthaler I (12004 ). assuming that 
the supermassive black hole is stationary at the Galactic center. These values for Uq and 
Wq are in good agreement with the Hipparcos results. 

The Hipparcos data used to determine Vq (the solar motion component in the direction 
of Galactic rotation) clearly show the well known "asymmetric drift," which when extrap- 
olated to zero dispersion should define the Local Standard of Rest (LSR). Our value of 
= 20 km s~^, based on HMSFR parallaxes and proper motions, is far above the asym- 
metric drift line, indicating that the HMSFRs as a group are orbiting the Galaxy slower 
than for circular orbits. Note that the youngest stars in the Hipparcos data, plotted at 
a dispersion of ~ 120 (km s~^)^, show a similar, but not as great a departure from the 
asymmetric drift line. Evidence that young stars lag the LSR orbit has also been found by 
Zabolotskikh. Rastorguev &: Dambis I (120021 ). 



Finally, we note that we find no evidence for a global motion of the LSR (i.e. disagree- 
ment with the Hipparcos solar motion) in the direction of the Galactic center or out of the 
plane of th e Galaxy larger than 6 k m s^^ (2cr). This is contrary to the conclusions of iKerr 



( 1l962l ) and lBlitz fc Spergel I (Il991al ). based on an analysis of HI data, that the LSR is moving 
away from the Galactic center at a speed of > 10 km s~^. 



3.2. Fundamental Galactic Parameters 

Since, as shown in §3.H HMSFRs are orbiting the Galaxy slower than for circular 
orbits, we must allow for such effects when modeling the Galaxy. In order to determine the 
fundamental parameters Rq and Bq, we solved for three additional parameters, allowing for 
an average peculiar motion for all sources with components Us toward the Galactic center 
(as seen by the source), Vg in the local direction of Galactic rotation and Wg toward the 
north Galactic pole. This solution, listed as Fit 2 in Table [3l yields Rq = 8.5 kpc and 
00 = 264 km s~^ and peculiar motion components of f/^, = 4 km s~^, Vg = —IQ km s~^ and 
Wg = 3 km s~^ The residuals show greatly reduced systematic deviations, and the value 
improved significantly to 112 for 67 degrees of freedom, compared to the solution without 
the average peculiar motions (Fit 1 in Table [3]). 

Two sources from the sample, NGC 7538 and G 23.6—0.1, displayed post-fit residuals 
significantly greater (> 3a) than the others. Removing these sources, we arrive at our 
"basic sample" of 16 HMSFRs. We repeated the fitting and found Rq = 8.40 ± 0.36 kpc. 
Go = 254 ± 16 km s-\ 17^ = 2.3 ± 2.1 km s-\ K = -14.7 ± 1.8 km s^^ and Wg = 
3.0 ± 2.1 km s^^ (see Fit 3 in Table [3]). The value for this sample was considerably 
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Fig. 4. — Solar motion components determined from Hipparcos stars (i.e. the reflex of the average 
motion of stars) versus stellar velocity dispersion after iDehnen &: Binney I (|l998l ). Top Panel: Vq is 
the Solar Motion in the direction of Galactic rotation (i.e. toward i = 90°). The "asymmetric drift" 
is shown with the dashed line. Middle Panel: Uq is the Solar Motion toward the Galactic center. 
Bottom Panel: Wq is toward the north Galactic pole. Also plotted at 50 (km s~^)'^ dispersion 
with open red squares are solar motion parameters obtained from the parallax and proper motions 
of star forming regions, and at zero dis persion with an open t riang le is the Wq component inferred 
from the proper motion of Sgr A* by iReid Brunthaler I (120041 ). Note the good agreement of 
the Uq and Wq components between Hipparcos and this study. The large deviation of the Vq 
component from the asymmetric drift from this study is not indicative of large V0 value, but points 
to a significant deviation from circular orbits for very young stars. 
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improved: 65.7 for 59 degrees of freedom. The near-zero average motion out of the plane 
of the Galaxy (Wg) is as expected for massive star forming regions. The residual motions 
in the plane of the Galaxy are shown in Figure O Most of the star forming regions have 
residual velocities consistent with measurement error combined with expected Virial motions 
within HMSFRs of ~ 7 km s~^ per coordinate. The most distant sources at low Declination 
(and low Galactic longitude) have larger residual velocities owing to greater parallax and 
proper motion measurement uncertainty and the scaling of proper motions to linear speeds 
by multiplying by distance. 

We feel that this solution provides the best estimates of the parameters for the current 
data set, under the assumption of a flat rotation curve. In §3.31 we show that the estimate 
of Rq is somewhat sensitive to the nature of the rotation curve of the Galaxy, leading to a 
systematic source of uncertainty for Rq of approximately ±0.5 kpc. Combining the statistical 
and systematic uncertainties in quadrature, we find Rq = 8.4 ± 0.6 kpc. 

The correlation coefficient between Rq and Go was 0.87, while all others were small. 
This is expected, since kinematic model distances increase with Rq and inversely with ©q. 
Thus, the ratio Qq/Rq, which is the angular rotation rate of the LSR, is determined to 
much better accuracy than either parameter separately. Holding Rq = 8.50 kpc (the lAU 
recommended value), we find 0o = 257.9 ± 7.7 km s~^ or Qq/Rq = 30.3 ± 0.9 km s~^ kpc^^. 
There is only a slight dependence of Qq/Rq on the value adopted for Rq. For example, setting 
Rq = 8.00 kpc, we obtain O0/-R0 = 30.0 ± 0.9 km s~^ kpc~^. See [JHlfor a discussion of the 
significance of this result. 

As shown in §3.3[ while estimates of Go change by ±20 km s~^ among the fits using 
different rotation curves, this variation can be accounted for mostly through the correlation 
with Rq, and, therefore, the least-squares fitting process incorporates this correlation in the 
formal uncertainty estimate. Thus, we conclude that the formal uncertainty of ±16 km s~^ 
for 00 is reasonable (provided that Rq is within 0.5 kpc of 8.4 kpc). When Rq is ultimately 
measured with much higher accuracy, Bq would be even better determined from the well 
determined ratio of 6o/-Ro- 

We also considered the possibility that a large positive value for Ug (toward the Galactic 
center), as could be expected from spiral density wave theory, might infiate the estimate of 
00. Holding Us = 17 km s^^(15 km s^^ greater than our best fit) did not significantly reduce 
the estimate of Oo, but did dramatically increase the to 200.1. Thus, we exclude a large 
Us value and that it could contribute to significant uncertainty in 0o. 
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Fig. 5. — Peculiar motion vectors of high mass star forming regions (superposed on an artist 
conception) after transforming to a reference frame rotating with the Galaxy, using best-fit values 
of Rq = 8.4 kpc and Qq = 254 km s"^ and removing an average motion of 15 km s~^ counter to 
Galactic rotation and 2 km toward the Galactic center. A 10 km s~^ motion scale is in the 
lower left. The Galaxy is viewed from the north Galactic pole and rotates clockwise. 



Table 3. Least-squares Fitting Results 



Fit 






eo 


dO 
dR 


Us 




Vs 












DF 






(kpc) 




(km s^^) 


(km s^^ kpc~^) 


(km s" 




(km s^^ 


) 


(km s^ 










(km s~^ kpc^^) 


1 


8.24±0 


55 


265±26 


0.0 


0.0 




0.0 




0.0 




263 


3 


70 


32.4±1.3 


2 


8.50±0 


44 


264±19 


0.0 


3.9±2 


5 


-15.9±2 


1 


3.1±2 


5 


111 


5 


67 


31.1±1.1 


3 


8.40±0 


36 


254±16 


0.0 


2.3±2 


1 


-14.7±1 


8 


3.0±2 


2 


66 


7 


59 


30.3±0.9 


4 


9.04±0 


44 


287±19 


2.3±0.9 


1.9±2 





-15.5±1 


7 


3.0±2 


1 


59 





58 


31.1±0.9 


5 


8.73±0 


37 


272±15 


Clemens- 10 


1.7±1 


9 


-12.2±1 


7 


3.1±1 


9 


52 


9 


59 


31.0±0.8 


6 


7.88±0 


30 


230±12 


Clemens-8.5 


2.7±2 


2 


-12.4±1 


9 


3.1±2 


3 


71 


2 


59 


29.6±1.0 


7 


8.79±0 


33 


275±13 


Brand-Blitz 


1.9±2 





-18.9±1 


8 


3.0±2 


1 


59 





59 


31.0±0.9 



Note. — Fits 1 & 2 used all 18 sources in Table [T] and have high values, owing to two outliers: NGC 7538 and G 23.6-0.1. 
Fit 3 excludes the two outliers and provides our basic result, under the assumption of a flat rotation curve. Fits 4-7 explore 
the effects of non-fiat rotation curves. "DF" is the degrees of freedom for the fit (i.e. number of data equations minus number of 
parameters). ( Us,Vs,Ws) are average pec uliar motions common to all sources (see Tableland Figure [8|), assuming the Hipparcos 
solar motion of Dehnen fc Binnev ( 1998h (see discussion in ij3.ll). All 6 n /i^n e stimates were obtained by holding Rq = 8.50 kpc and 

rotation curves for (i?o[kpc],0o[km s~^]) = (10,250) 
rotation curve. Both the Clemens and Brand-Blitz 



solving for Qq- "Clemens-10" and "Clemens-8.5" refer t o the [ Clemens 
and (8.5,220), respectively; "Brand-Blitz" refers to the iBrand &: Blitz 




rotation curves were scaled to the fitted values of Rq and Oq. 
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3.3. Rotation Curves 



We have until now assumed tliat tlie rotation curve of the Galaxy is flat (i.e. Q{R) = Go). 
In order to investigate deviations from a flat rotation curve, we used the basic sample of 



16 sources and added the parameter ^ to the model. 



A least-squares fit yielded ^ 



2.3 ± 0.9 km s~^ kpc^^, with an improved compared to the flat rotation curve fit (see 
Fit 4 in Table [3]), but with an increased correlation coefficient between Rq and ©o of 0.90. 
We tested how sensitive ^ was to the two outer Galaxy sources by dropping S 269 and 
WB 89-437 from the sample and re-fitting. This yielded ^ = 1.9 ± 1.2 km s"^ kpc~^ and 
indicated that these sources do not provide all the leverage for a rising rotation curve. Thus, 
we find a nearly flat rotation curve between Galactocentric radii of about 4 to 13 kpc, with 



conclusions reached in a number of papers ( 


Fich. Blitz & Stark 


1989; 


Brand & Blitz 


1993; 


Honma & Sofue 


1997; 


Maciel & Laeo 


2005h. For example 


Fich. Blitz & Stark 


( 


1989 


) find 



that the rotation curve is nearly flat for Gq 
Gn = 250 km s-\ 



220 km s ^ and that it rises gradually for 



We also tested more complex rotation curves b y re placing the simp l e linea r form just dis- 
cussed with the rotation curves of lClemens I (119851 ) and lBrand &: Blitz I ( 1l993l ). Clemens sup- 
plied two curves: one assuming the old lAU constants of Rq = 10 kpc and Go = 250 km s~^ 
and the other assuming the new constants of -Ro = 8-5 kpc and Go = 220 km s^^. These 
models have slightly different shapes, with the old model generally having rotational speeds 
that rise faster with radius than the new model. For either model, we fitted for different 
values of Rq (which we used to scale model radii) and Go (which we used to scale rotation 
speeds). The fit using the old model, listed as Fit 5 in Table [31 gave Rq = 8.7 ± 0.4 kpc 
and Go = 272 ± 15 km s~^, with an improved = 52.9 for 59 degrees of freedom compared 
to our solution for a flat rotation curve. The improvement is partly from a better match to 
the two sources in the Outer arm (S 269 and WB 89-437). Using the new rotation model, 
gave Rq = 7.9 ± 0.3 kpc and Go = 230 ± 12 km s ~^, w ith a considerably worse = 71.2 
(see Fit 6 in Table [3]). Using the iBrand fc Blitz I (1l993l ) rotation curve, also scaled by the 
fitted values of Rq and Go, we obtain Fit 7 in Tabled with values of i?o = 8.8 ± 0.4 kpc and 
Go = 275 ± 15 km s~^ and a = 59.0 for 59 degrees of freedom, intermediate between the 
values for the two Clemens models. 

Clearly there is some sensitivity of the best fit -Ro value to the models, and we adopt a 
systematic uncertainty in Rq of ±0.5 kpc. Note that, as discussed in §3.21 the ratio G0/-R0 
has much less modeling sensitivity. With the current parallax and proper motion data, 
we cannot conclusively distinguish among the rotation curves presented. However, with the 
many more parallaxes and proper motions expected in the next few years from the VLBA and 
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VERA telescopes, we should be able to make considerable progress in refining the rotation 
curve of the Milky Way. 

4. Kinematic Distances 

Figure E] compares the locations of the star forming regions determined by trigonometric 
parallax and by kinematic distances. The kinematic distances were computed for the lAU 
standards Rq = 8.5 kpc and Oq = 220 km s~^ and the standard definition of LSR. For 13 
of 18 regions (11 of 16 in the basic sample), the kinematic distance exceeds the true source 
distance; in 3 cases the discrepancy is over a factor of two. The kinematic distances for 
(these) star forming regions tend to over-estimate the source distances. 

As shown above, HMSFRs on average orbit the Galaxy ^15 km s~^ slower than the 
circular rotation speed. Taking this into account, a prescription for a "revised" kinematic 
distance for a high-mass star forming region could be as follows: 

1) add back the (old) Standard Solar Motion corrections to the LSR velocities, returning 
them to the heliocentric frame; 

2) apply "best values" for the solar motion to calculate a revised "LSR" velocity, v^^g-^; 

3) subtract —15 km s^^ from the velocity component in the direction of Galactic rotation; 

4) calculate a kinematic distance using values for the fundamental parameters of the Milky 
Way, e.g. Ro = 8.4 kpc and Go = 254 km s~^, that are consistent with astrometric 
measurements; and 

5) when determining the uncertainty in the kinematic distance, include a systematic con- 
tribution allowing for the possibility of a 7 km s~^ uncertainty in f^g^. 

Table H] gives parallax distances, standard (old) kinematic distances, and revised kine- 
matic distances and uncertainties (using the above prescription) for all 18 HMSFRs listed 
in Table [H (We provide the FORTRAN source code used to calculate revised kinematic 
distances in the on-line material.) Note that our prescription for the uncertainty in kine- 
matic distances performs reasonably well for our basic sample (excluding the two sources 
G 23.6—0.1 and NGC 7538 which we earlier noted as outliers). The mean difference between 
the parallax and kinematic distances is near zero and the differences divided by their uncer- 
tainties average to near unity. Now only half (8 of 16) of the sources in the basic sample 
have kinematic distances that exceed the true source distance. There are no cases for which 
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Fig. 6. — The locations of the star forming regions determined by trigonometric parahax {dark 
blue circles) and by kinematic distances {light magenta circles), assuming lAU recommended values 
of i?o = 8.5 kpc and 0o = 220 km s~^ and the Standard Solar Motion to define the LSR. 
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the discrepancy is a factor of two, and the estimated uncertainties reasonably account for 
differences between the parallax and kinematic distances. 

While the prescription outlined above results in some improvement in kinematic dis- 
tances compared to the standard approach, the improvement is not as great as one might at 
first expect. This occurs because the definition of the LSR uses the Standard Solar Motion. 
While the Standard Solar Motion differs only slightly from the Hipparcos solar motion for 
components toward the Galactic center (Uq) and the north Galactic pole (Wq), there is a 
large discrepancy for the component in the direction of Galactic rotation. The Standard 
value is Vq^'^ = 15.3 km s~^, whereas the Hipparcos value is Vq = 5.25 km s~^. The 
+10 km s~^ "error" in V^*'^ partially compensates for the 15 km s~^ slower Galactic orbits 
of HMSFRs shown in §3.1[ (Note that a positive change in the solar motion component 
Vq results in a negative change in a source peculiar motion component Vg). Even with 
the improved prescription for kinematic distances, one cannot really hope to discern spiral 
structure using kinematic distances. 



Galactic Coordinates 



There is excellent agreement between the two independent VLBI measurements of 
Oo/-Ro- the measurement based on parallaxes a nd proper motions of HM SFRs (this pa- 
per) and based on the proper motion of Sgr A* (IReid fc Brunthaler II2004I ). This gives us 
confidence that 1) we can well model the Galaxy with parallax and proper motions of HMS- 
FRs and 2) Sgr A* is indeed a supermassive black hole at the dynamical center of the Milky 
Way. These findings offer an independent definition of the Galactic plane and Galactic co- 



actic plane is based primarily on the thin 
1960). The Sun is defined 



ordinates. Currently, the lAU definition of the Gal 
distribution of neutral hydrogen 21 cm emission (IBlaauw et al. 
to be precisely in the plane and the origin of longitude was set by the centroid of the radio 
emission of th e large , comp lex source Sgr A. The Sun is now known to be ~ 20 pc north of 
the plane (see iReed I (120061 ) and references therein) and the supermassive black hole, Sgr A*, 
is offset by a few arcmin from the lAU defined center. 

In the future, one could consider redefining Galactic coordinates based, in part, on 
the proper motion of Sgr A*, which, after correction for the well- determined solar motion 
component perpendicular to the Galactic plane, gives the orbital plane of the Local Standard 
of Rest. The zero of longitude would be best defined by the position of Sgr A*. This would 
place our supermassive black hole at the origin of Galactic coordinates, and one could rotate 
the reference frame to remove the Sun from its special location precisely in the Galactic 
plane. 
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Table 4. Parallaxes vs. Kinematic Distances 



Source 






















(Ac 


rr\ 
gj 




V) 


km s ^ 


(kpc) 


(kpc) 


(kpc) 


G 23.0-0.4... 


23. 


,01 


-0. 


,41 


+81 


4.59 


4.97 


4.72+H 

— U.ci 


G 23.4-0.2... 


23. 


,44 


-0, 


.18 


+97 


5.88 


5.60 


5.29+°-^ 

^ ^ —0.6 


G 23.6-0.1... 


23. 


,66 


-0, 


.13 


+83 


3.19 


5.04 


4.77+n? 

— U.cS 


G 35.2-0.7... 


35. 


,20 


-0, 


.74 


+28 


2.19 


2.00 


1.99+°-^ 


G 35.2-1.7... 


35. 


,20 


-1, 


.74 


+42 


3.27 


2.85 


2.76+°-! 

—0.4 


W 51 IRS 2.. 


49. 


,49 


-0, 


.37 


+56 


5.13 


5.52 




G 59.7+0.1... 


59. 


,78 


+0, 


.06 


+27 


2.16 


3.07 


3.45+il 


Cep A 


109. 


,87 


+2, 


.11 


-10 


0.70 


1.09 


0.551°:^^ 


NGC 7538.... 


111. 


,54 


+0, 


.78 


-57 


2.65 


5.61 




IRAS 00420.. 


122. 


,02 


-7, 


.07 


-44 


2.13 


3.97 


3.181^ 


NGC 281 


123. 


,07 


-6 


.27 


-31 


2.82 


2.69 


2.081°:^ 


W3(0H) 


133. 


,95 


+1, 


.06 


-45 


1.95 


4.28 


3.421°:^ 


WB 89-437... 


135. 


,28 


+2, 


.80 


-72 


5.99 


8.68 


6.8911:^ 


S 252 


188. 


,95 


+0, 


.89 


+ 11 


2.10 


4.06 


o oq+4.2 
O . OO 2 4 


S 269 


196. 


,45 


-1, 


.68 


+20 


5.29 


4.13 


3.351?:° 


Orion 


209. 


,01 


-19 


.38 


+10 


0.41 


0.99 


71 +°-'^ 


G 232.6+1.0.. 


232. 


,62 


+ 1. 


.00 


+23 


1.68 


1.92 


1.44l°| 


VY CMa 


239. 


,35 


-5, 


.06 


+18 


1.14 


1.56 


1 10+°-^ 



Note. — is the measured parallax converted to distance; -Of is the 
kinematic distance based on standard LSR velocities; D^"^" and a{Dk) is the 
revised kinematic distance and its uncertainty, calculated for Rq = 8.4 kpc, Oq = 
254 km s~^ and Us = —15 km s^^, following the prescription outlined in fJH All 
kinematic distances assume a flat rotation curve. 



- 25 - 



Discussion 



Very Long Baseline Interferometry now routinely yields parallax measurements with 
accuracies of ~ 10 /xas, corresponding to 10% uncertainty at a distance of 10 kpc, and proper 
motions that are usually accurate to ~ 10 /xas y~^ or better than than ~ 1 km s~^ at similar 
distances. Target sources include molecular masers associated with star formation and red 
giant stars, as well as non-thermal continuum emission associated with young T Tau stars 
and cool dwarfs. Combining the first results of parallaxes for high-mass star forming regions 
from the VLB A and the Japanese VERA project has allowed us to begin to investigate the 
spiral structure and kinematics of the Galaxy. 

We have accurately located three of the spiral arms of the Milky Way and directly mea- 
sured a pitch angle of 16° for a portion of the Perseus sp iral arm. This pit ch angle is similar to 
those of spiral arms in other galaxies of type Sb to Sc (IKennicutt Ill98ll ). Two armed spirals 
can account for most of the known large H II regions only if the arms wrap twice around the 
Galaxy; this requires pitch angles of ~ 8°. With a pitch angles greater than ^12°, the Galaxy 
needs to have four arms in order to account for the approximate locations of H II regions 
( iGeorgelin fc Georgelin Ill976l : iTaylor fc Cordes Ill993l ). There has been co nsiderable discus- 



sion in the l iterature concerning the number of spiral arms in the Galaxv (ISimonson 



Bash 



1981 



Vallee 



1995 



Drimmel 



2000 



Drimmel fc Spergel 



2001 



Benjamin et al. 



1976 



2005 



Nakanishi fc Sofue 



2006; Steiman-Camerson. Wolfire fc HoUenbach 



20081), with Spitzer GLIMPSE 



survey results sug gesting only two arms can be traced in the redder, older population of stars 
( iBenjamin I l2008l ). Perhaps the VLBI and infrared survey results can be reconciled if the 
Milky Way exhibits a hybrid structure, consisting of two dominant spiral arms, populated 
by both young and old stars and with pitch angles near 16°, and two weaker arms traced 
only by young stars. 

Our finding that HMSFRs on average orbit the Galaxy 15 km s~^ slower than expected 
for circular orbits has implications for star formation and spiral density wave theory. The 
plot of the apparent solar motion in the direction of Galactic rotation (Vq) versus stellar 
dispersion (Figure Hj) can be interpreted as a time sequence, with stellar age increasing with 
dispersion. The 15 km s~^ slower orbital speed of HMSFRs displays as a positive departure 
of the apparent Solar Motion with respect to the asymmetric drift (the fitted trend to the 
Hipparcos data shown in Figure H]), since the Sun appears to orbit faster when measured 
against such stars. One explanation for this finding is that HMSFRs are born in elliptical 
Galactic orbits, near apocenter, with orbital eccentricity of about 0.06. As young stars 
continue to orbit the Galaxy, their orbits become more circularized, as evidenced by the lesser 
departure of the youngest Hipparcos bin (mostly late B-type stars) from the asymmetric drift 
line compared to the HMSFRs. The gradual transfer of angular momentum from gas to stars 
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in the Galaxy proposed by IChakrabarti I (120091 ) may explain this. At a stellar dispersion of 
~ 300 (km s~^)^, which corresponds to A2- to A5-type stars with colors B — V = 0.1 and 
characteristic main-sequence lifetimes of ~ 1 Gy, the stars join the asymmetric drift. As 
stars continue to age, their orbits are progressively "randomized" and they (again) become 
part of a slower orbiting population, which appears as a larger apparent Vq. 

Parallaxes measurements alone generally cannot yield Rq (except for a parallax of 
Sgr A*). However, since galaxies rotate in a fairly smooth fashion, a kinematic model 
can be directly compared with distance and relative motion measurements in order to es- 
timate Rq and 00. In this paper, we have demonstrated that parallax and proper motion 
measurements for HMSFRs across large portions of the Galaxy can separate estimates of Rq 
and 9o, although because of the somewhat restricted coverage of the Galaxy currently avail- 
able, we have a significant correlation between these parameters. Our best estimate of Rq is 
8.4 ± 0.36 ± 0.5 kpc, where the second uncertainty is systematic and comes from our lack of 
detailed knowledge of the rotation curve of the Galaxy. This estimate is consi stent with th e 
"best" Rq of 8.0 ± 0.5 kpc from a combination of many methods reviewed by iReid I (119931 ). 
Also, recent direct estimates of Rq from radial velocitie s and elliptical pa ths of stars that 
orbit Sgr A* have cori verged on values of 8.4 ± 0.4 kpc (iGhez et al. II2008I ) and 8.33 ± 0.35 
(iGillessen et al. 1 120091 ). (These estimates assume that Sgr A* is nearly motionless at the 
Galactic center. Relaxing this assumption decreases the estimates to about 8.0 kpc.) Of 
course, many other less direct estimates of Rq can be found in the literature and span a 
much greater range. 

The characteristic rotation speed of the Galaxy (Oq) is a crucial parameter not only for 
Galactic dynamics and kinematic distance determinations, but also for estimati ng the total 



mass in dark matter arid the history and fate of the Local Group of galaxies (ILoeb et al. 



20051 : IShattow fc Loeb I l2008l ). Estimates of the rotation spee d of the Galaxy from the 
recent literature span a very large ra nge between 184 km s~^ (jOUing fc Merrifield I Il998l ) 



and 272 km s~^ (iMendez et al. Ill999l ). Most estimates of Qq are based on analyses of the 
shear and vorticity of large samples of stars in the (extended) solar neighborhood and thus 
really measure Oort's A and B parameters. Quoted values of Qq then come by assuming a 
value for Rq and using the relation Qq = Rq{A — B) . Our result that Qq = 254 ± 16 km s~^ 
was obtained by fitting for both Rq and Qq using full 3-dimensional locations and motions 
of sources well beyond the extended solar neighborhood and, thus, does not assume a value 
for Rq. However, as discussed in §3.21 with the present distribution of sources there is 
considerable correlation between Rq and Qq parameters, which is reflected in the ±16 km s~^ 
formal uncertainty for Qq. 



Our estimate of the ratio 0o/-Ro of 30.3 ± 0.9 km s ^ kpc ^ is determined more accu- 
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rately than either parameter individually and is nearly independent of the value of Rq over 
the range of likely values between about 8.0 and 8.5 kpc. This value differs considerably from 
that determined from the lAU yalues of 0n/-Rn = 220 km s~^/8.5 kpc = 25.9 km s~^ kpc~^ 



and differs marginally from the iFeast fc Whitelock I (119971 ) analysis of Hipparcos Cepheids 
of 27. 19 ±0.87 km s~^ kpc~^. Recent studies, using samples of OB-type stars within 3 kpc of 



metrically determined distan ces), arrive at A 
2000h and 32 km kpc^^ 



B between 30 km s ^ kpc ^ ( 


Uemura et al. 


1998: 


Elias. Alfaro & Cabrera-Cano 


2006h. 



with uncertainties of about ±1.5 km s ^. 

Our value for G0/-R0 is in excellent agreement with that determined directly from the 
apparent proper motion of Sgr A* (the superrnassiy e black hole at the Galactic center) 
of 6.379 ± 0.024 mas y~^ (IReid fc Brunthaler II2004I ). One expects a supermassive black 
hole to be stationary at the dynamical center of the Galaxy to better than ~ 1 km s~^ 



(IChatterjee. Hernquist fc Loeb II2OO2I : iDorband. Hemsendorf fc Merritt II2003I : iReid fc Brunthaler 
20041 ). Hence, Sgr A*'s apparent motion should be dominated by the effects of the Galactic 
orbit of the Sun. After correcting for the Solar Motion of 5.25 km s~^ in the direction of 
Galactic rotation (jPehnen fc Binney 1 119981 ) . Sgr A*'s apparent motion yields a global esti- 
mate of Qq/Rq = 29.45 ± 0.15 km s~^ kpc~^. Thus, there is excellent agreement between 
this and our global and direct method for measuring Qq/Rq. 

Coupling the Sgr A* motion result of Qq/Rq = 29.45±0.15 km kpc~^ with estimates 
of -Ro from stellar orbits in the Galactic center (IGhez et al. I l2008l : iGillessen et al. I l2009l ) 
of 8.4 ± 0.4 kpc yields Go = 247 ± 12 km s~^. This result is also a direct and global 
measurement of Oq and is independent of our result from parallaxes and proper motions of 
star forming regions. Combining the Galactic center and star forming region estimates gives 
Go = 250 ± 10 km s~\ 

It seems clear that Go is near the upper end of the range of estimates in the literature. We 
note that both the Galactic center stellar orbit and the star forming region parallax results 
assume the Hipparcos solar motion of 5.25 km s~^ in the direction of Galactic rotation. Only 
if the interpretation of the asymmetric drift is incorrect or if the entire solar neighborhood 
orbits the Galactic center ~ 30 km s~^ slower than the Galaxy spins could Gq be equal to 
the lAU recommended value of 220 km s~^. 

We have determined the rotation speed of the Milky Way at the radius of the Sun to be 
^ 250 km s~^ and the rotation curve to be nearly flat or slightly rising with distance from the 
Galactic center. These values are nearly identical to those of the Andromeda galaxy (M31) 
as shown in Fig u re [71 The rotation curve of Andromeda, deterr nined from HI emission by 
Carignan et al. I (120061 ) based on interferometric observations of lUnwinI (119831 ) . indicates a 



- 28 - 



300 



250 



200 



150 



100 



Milky Way 




□ □ 



Andromeda (M 31^ 



10 



15 



20 



R (kpc) 



Fig. 7. — Rotation speed versus radius for the Andromeda gal axy and the Milky Wa y. The red 



squares are based on HI observations of Andromeda tabulated by lCarignan et al. I (|2006l ). The blue 
filled circle is our best estimate of Gq = 254 it 16 km s~^ at Rq = 8.4 kpc for the Milky Way, derived 
from the parallax and proper motions of high mass star forming regions. The blue dot-dashed line 
is for a flat rotation curve, and the blue dashed line corresponds to a slightly rising rotation curve of 
2.3 km s~^ kpc~^(see ^3.3p . These lines are plotted over the range of Galactocentric radii sampled 
by the parallax and proper motion results. Note that these two galaxies have nearly identical 
rotation speeds over this range. 
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speed of 251 km at a radius of 8 kpc, a slightly rising curve out to about 15 kpc, and a 
slow dropoff to about 225 km s~^ beyond 20 kpc. The most straightforward interpretation 
of the similarities of the rotation curves for the Milky Way and Andromeda is that these two 
galaxies are nearly equal in size and mass. 



Finally, we note that iReid fc Brunthaler I (120041 ) placed a strong upper limit of —0.4 ± 
0.9 km s~^ for the component of peculiar motion of Sgr A* perpendicular to the plane of the 
Galaxy. However, the determination of the component in the direction of Galactic rotation 
was considerably less accurate: 18 ± 7 km s~^, as one must remove the uncertain effects of 
the solar orbit. Reid & Brunthaler did this by removing 27 .19±0.87 km s~^ kp c ~^, ba sed on 
Hipparcos measurements of Oort's constants {A — B) by iFeast fc Whitelock I (119971 ). from 
the observed motion of Sgr A* in the Galactic plane of 29.45 km s~^ kpc~^. This method 
assumes that Qq/ Rq = A — B and that estimates of the shear and vorticity of nearby stars 
from Hipparcos data indicate the large-scale differential rotation of the Galaxy and are not 
subject to local irregularities in the solar neighborhood. Since we now have a direct, global 
estimate of 0o/-Ro = 30.3 ± 0.9 km s~^ kpc~^, we find the peculiar motion of Sgr A* in the 
direction of Galactic rotation to be —7.2 ± 8.5 km s~^, with little sensitivity to Rq (adopted 
to be 8.5 kpc here). This adds additional strong evidence that Sgr A* is a supermassive 
black hole, which is nearly stationary at the dynamical center of the Galaxy. 
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7. Appendix 

Since we have measured the position, distance, LSR velocity and proper motion of each 
source, we know its full 3-dimensional location in the Galaxy and full space motion relative 
to the Sun. Given a model of Galactic rotation, we can then calculate the non-circular 
(peculiar) velocity of each source. While this calculation is conceptually simple, in practice, 
there are some subtleties and sign convention issues that can lead to errors, and so here we 
present the necessary formulae (and FORTRAN source code in the on-line material). 

The required Galactic and Solar Motion parameters are given in Table and those 
associated with the source are defined in Table O A schematic depiction of these parameters 
is given in Figure [H We assume that the Sun is in the Galactic plane and calculate a source's 
peculiar motion (i.e. with respect to a circular Galactic orbit) as follows: 

We convert wlsr to a heliocentric frame, vneiio, by adding back the component of the 
Standard Solar Motion in the line-of-sight direction that had been removed from the ob- 
served Doppler shift to calculate wlsr- Note that one needs to use the (old) Standard Solar 
Motion, which defines the LSR frame, and not the best values available today. Generally, ob- 
servatories have adopted a value of 20 km s"^ toward a(1900) = 18*^, 5(1900) = +30*^ for the 
Standard Solar Motion. Precessing these coordinates to the epoch of observation 2006) 
and converting to Galactic Cartesian coordinates yields the {Uq'^,Vq^'^,Wq''^'^) values listed 
in Table [51 Then, 

VHeiio = vlsr - (^l*'^ COS i + V^^'^ siu i) COS b - W^^'^ sin b . 

Rotate the motion vector from the equatorial heliocentric frame {fJ'a,fi'S,VHeiio) to a 
Galactic heliocentric frame [fiijfibjVHeHo)- This is a rota.tion a bout a radial axis and is de- 



fined by the lAU in B195 coordinates (IBlaauw et al. I Il960l ). For coordinates in J2000 



Reid fc Brunthaler I (12004 ) give the right ascension and declination of the NGP as ap = 
12'^51™26f2817 and 6p = 27°07'42'.'013, respectively, and the zero of longitude is the great 
semicircle originating at the NGP at the position angle 6 = 122?932. Galactic latitude can 
be obtained from 

sin 6 = sin 5 cos (90° — 5p) — cos 5 sin {a — ap — 6^*) sin (90° — Sp) . 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



Table 5. Galactic and Solar Parameters and Nominal Values 



Parameter Value Definition 

Rq 8.5 kpc Distance to the GC (lAU value) 

00 220 km s^^ Rotation speed of LSR (lAU value) 

Qs 220 km s~^ Rotation speed of Galaxy at source 



10.3 km s"^ Standard Solar Motion toward GC 

V^^"^ 15.3 km s'^ Standard Solar Motion toward £ = 90° 

7.7 km s-^ Standard Solar Motion toward NGP 

Uq 10.0 km s~^ Hipparcos Solar Motion toward GC 

Vq 5.2 km s~^ Hipparcos Solar Motion toward i = 90° 

Wq 7.2 km s^^ Hipparcos Solar Motion toward NGP 



Note. — GC: the Galactic Center; LSR: Local Standard of Rest; NGP: 
North Galactic Pole. The Standard Solar Motion must be used to convert from 
VhSR to VHeiio, sincc (hopefully) all observatories have used this definition. The 
values given above come from an assumed Solar Motion of 20 km s~^ toward 
R.A.(1900)=18^ an d Dec.(1 900)-30° processed to J2000.0. Hipparcos Solar Mo- 
tion values are from lDehnen fc Binnev i (|1998i) . 
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Table 6. Source Parameter Definitions 



Parameter Definition 



£ Galactic longitude 

D Distance from Sun (I/tTs) 

Dp Distance from Sun projected in plane 

Rp Distance from GC projected in plane 

wlsr ■■■■ LSR radial velocity 

VHeiio ■■■ Heliocentric radial velocity 

/Xq, Proper motion in R.A. (/i^ = fix/ cos 5) 

l^s Proper motion in Decl. (/x^ = fj.y) 

P Angle: Sun-GC-source 

Us Peculiar motion locally toward GC 

Vs Peculiar motion locally in direction of Galactic rotation 

W. Peculiar motion toward NGP 



Note. — GC: the Galactic Center; LSR: Local Standard of Rest; NGP: North 
Galactic Pole. 
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Fig. 8. — Schematic depiction of source and Galactic parameters. 

A useful angle (f) can be determined (between 0°and 360°) from 

smcj) = ^cos5sin {a — ap — 6'*) cos (90° — 5p) + sin 5 sin (90° — 5p)^ / cos 6 

and 

cos = cos 5 cos {a — ap — Q'^) / cos b , 
and then Galactic longitude follows from 

£ = 0+(^-9O°) . 

Proper motion in Galactic coordinates {fii,fiiy) can be easily calculated from the motion 
in equatorial coordinates (/Xq, fis) by differencing (£, b) values for coordinates determined, 
say, one year apart. This usually requires 64-bit precision in the calculations. Note that /i/ 
will naturally be defined positive in the direction of increasing Galactic longitude, which is 
counter to Galactic rotation. 



Convert the proper motions to linear speeds (by multiplying by distance) via 



ve = Djii cos 6 and Vh = Dji}, , 
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where ni cos b is the actual motion tangent to the sky in the direction of Galactic longitude. 
We now convert from spherical to Cartesian Galactic coordinates at the location of the 

Sun. 

Ui — {vHeiio COS b — Vb sin b) cos i — sin £ , 

Vi = {vHeiio COS b — Vb sin b) sin i + V£ cos i , 
Wi = Vb cos b + VHeiio sin b . 

Next add the full orbital motion of the Sun, using the best values of the solar motion 
and the circular rotation of the Galaxy at the position of the Sun (U^ yQ+Qo,WQ): 

U2^Ui + u^ , V2 = Vi + v^ + eo , W2^Wi + w^ . 

The Galactocentric distance to the source projected onto the Galactic plane is given by 

R/ ^ Ro'^ + Dl-2RoDpCOse , 

where Dp = D cos b. The angle f3 between the Sun and the source as viewed from the 
Galactic center can be determined in all cases (i.e. from 0° to 360°) from 

sin (3 = — ^ sin i and cos P — — — p^^^^ ^ 

Rotate the vector {U2, V2, W2) through the angle P in the plane of the Galaxy and remove 
circular Galactic rotation at the location of the source to yield (Us, Vg, Wg) : 

Us — U2 cos /9 — V2 sin /9 , 

K = y2Cos/3 + C/2sin/3 -0, , 

The vector {Us^Vg^Ws) gives the non-circular (peculiar) motion of the source in a Carte- 
sian Galactocentric frame, where Ug is radially inward toward the Galactic center (as viewed 
by the source), Vg is in the local direction of Galactic rotation and Wg is toward the north 
Galactic pole. 



